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The paper presents a certainty equivalence output feedback backstepping adaptive con- 
trol design method for the systems of any relative degree with unmatched uncertainties 
without over-parametrization. It uses a fast prediction model to estimate the unknown 
parameters, which is independent of the control design. It is shown that the system’s in- 
put and output tracking errors can be systematically decreased by the proper choice of the 
design parameters. The approach is applied to aerospace control problems and tested in 
numerical simulations. 


I. Introduction 

Adaptive control problems are challenging for systems with unmatched uncertainties, and the majority of 
design methods is based on the backstepping technique outlined in Ref . 6 They become more complex when 
the whole state of the system is not available for feedback. It has been shown in Ref . 6 that direct application 
of the certainty equivalence principle leads to over-parametrization. To avoid it, nonlinear damping terms 
were introduced in the control design, which leads to the adaptation rate to enter into the control law. This 
may result in the high magnitude control signals in the case of large adaptive rates (fast adaptation) , which is 
desirable from the perspective of the estimation of unknown parameters. An alternative certainty equivalence 
control design method for state feedback problems is presented in Ref ., * 1 which avoids over parametrization 
for linear systems, but the method is only applicable for nonlinear systems with the relative degree not 
exceeding two. 

In Ref ., 11 we introduced a certainty equivalence state feedback indirect adaptive control approach without 
over parametrization for nonlinear systems of any relative degree. The approach was based on the state 
prediction model, which is capable of providing fast estimation of unknown parameters independent of the 
control design. This property is the consequence of feeding back an error term with the gain proportional to 
the square root of the adaptation rate, like in the modified reference model MRAC (M-MRAC) architecture 
introduced in Ref ., 10 where it has been shown that the error feedback gain acts as a damping factor for the 
adaptive signals, whereas the adaptation rate determines its frequency. In this paper we extend the method 
to the systems in parametric output feedback form. This extension uses conventional filtering approach 
to transform the system into the form suitable for the parameter estimation. The control design uses the 
prediction model’s and the filter’s states, and follows the command filtered backstepping procedure from 
Ref . 2 It is shown that the input tracking error (difference between ideal control and command filtered 
certainty equivalence control signal) and output tracking error can be regulated as desired by the proper 
choice of design parameters. 

The rest of the paper is organized as follows. In Section II, we give the problem statement, main 
assumptions and input-output filtered transformations. In Section III, we introduce the identification model 
and give its properties. The control design and performance analysis are presented in Section IV. Aerospace 
applications are presented in Section V, and some concluding remarks are given in Section VI. 
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II. Problem Statement 


Consider an uncertain single input single output (SISO) system in the parametric output feedback form 6 
(p. 99) 


x(t) = A n x(t) + f(y) + $(y)6 + du(t) 

y(t) = cZx{t), (i) 

with some initial conditions a:(0) = Xq, where x € R n , u £ R and y € R are the state, input and output 
of the system, 9 £ R p is a vector of unknown constant parameters, f : R —> R" and $ : R — 1 R nxp are 
sufficiently smooth known functions, and A n , d , c n have forms 
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It is assumed that d is unknown, d m is positive and bounded below by a known constant e?*, and the 
polynomial d m s m + ■ ■ ■ + d\s + do is known to be Hurwitz. 

The objective is to design a control signal u(t) using only the output y(t) such that all closed-loop signals 
are bounded, and y(t) tracks the output y r (t) = c T x r (t) of the reference model 

x r (t) = Ax r (t) + br(t), a; r (0) = x 0 , (2) 

where x r £ R 9 (q = n — m) is the state of the model, r(t) is a piece-wise continuous and bounded external 
command, A = A q — b q k T , b = kb q . Here we denote b q = [0( g _!) Xl 1] T , and the gains k and k are chosen 
to make A Hurwitz and meet the performance specifications. 

Following the steps from Ref. 6 (p.329), we transform the system (1) by means of the filters introduced 
in Ref. 5 In this particular case these filters take the form 

s(t) = A 0 s(t) + k 0 y(t) + f(y) 

E(t) = A 0 E(t) + $(y) 

l±(t) = Aofjb(t) + b n u(t) , (3) 

where the constant vector fco is chosen such that the matrix Ho = A n — k^c^ is Hurwitz. It is straightforward 
to show that 

m 

x(t) = s(t) + E(t)0 + ^ > ( 4 ) 

2—0 

where Wi(t) = Hg/i.(t), * = 0, . . . ,m, and S(t) is a solution of the exponentially stable system 

m = AoS(t ) . (5) 

In Ref., 6 5(t) was treated as a bounded external disturbance, and a nonlinear damping term was added in 
the control signal to compensate for it. Here, we show that the certainty equivalent controller can achieve 
the objective without extra terms or over-parametrization. 

Using the representation (4) we transform the output dynamics as follows 

y(t) = x 2 (t) + fAy) + $(i)( 2/)0 

m 

= s 2 (t) + 3(2 )(t)0 + ^2 diWi i2 (t) + S 2 (t ) + fi(y) + 

*= o 

= d m w m , 2 (t) + fi(y) + s 2 (t) + h T (t)d + 6 2 (t ) , (6) 

where S ( 2 )(t) denotes the second column of matrix 3(f), = [d m _i ... do 0 T ] T is the augmented vector 

of unknown parameters, and h(t) = [w m -i, 2 (t) ■ ■ ■ Wo <2 {t ) S( 2 )(f) + $( 1 )(y)] T is the augmented vector of 
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regressors. Combining (6) with the iD m -dynamics we obtain a g-dimensional system for the control design 

y(t) = d m Wm,2(t) + fi{y) + s 2 (t) + h T (t)tf + S 2 (t) 

(f) — 'U- ! rn.j+l (^) kjW m ,l if ) : J — 2, . . . , Q 1 

Wm,q{t ) = U(t) + W m ,q+l(t) - k q W m ,l(t) . (7) 

Obviously (7) is in parametric output feedback form, the states of which are available for the control de- 
sign. Therefore, the state feedback M-MRAC backstepping approach from the Ref. 11 can be applied. The 
difference is the presence of the unknown virtual control confficient d m . Since the control law will contain 
division by the estimate of d m , a projection operator can be employed to bound it away from zero using the 
available lower limit d*. 


III. Identification Model 

Since the uncertainties are lumped in the first equation of (7), we use a reduced order identification model 
to estimate the unknown parameters 

y{t) = dmit)w m ^{t) + fi(y) + s 2 (t) + /i T (t)i9(t) + A y{t) 

2/(0) = 2/o , (8) 

where y{t) is the output prediction, y(t ) = y(t) — y(t) is the output prediction error, A > 0 is a design 
parameter, d m (t) and $(t) are the estimates of the unknown parameters, which are generated according to 
adaptive laws 

dmit) = 7 Pr (d m (f),y(t)w m , 2 (0) (9) 

= iy(t)h(t ) , 

where 7 > 0 is the adaptation rate and Pr (•, •) denotes the projection operator (see Ref. 8 for details). 

The output prediction error dynamics do not explicitly depend on the control signal 

y(t) = -Xy(t) + rj(t) + 6 2 (t), (10) 

where we define fj(t) = y(t) - fj(t) = d m w m , 2 (t) + h T (t)tf - dm(t)w mi2 (t ) - /r T (t)i9(t) = dm{t)w mi2 (t) + 
h 1 {t)'d{t) with d m (t) = d m — d m (t) and i)(f) = i? — 1 9(t) being the parameter estimation errors. 

The following properties of the identification model can be established. 

Lemma III.l The error signals y{t), d m (t ) and id(t) are globally uniformly bounded, and y(t) — > 0 as 
t — > 00 . 

Proof. Consider a candidate Lyapunov function 

V(t ) = y 2 (t) + -d 2 m (t) + V(t)tf(t) + 6 r (t)P6(t) , (11) 

7 7 

where P is the positive definite solution of the Lyapunov equation A J P + PA 0 = —Q for a given positive 
definite matrix Q. It is straightforward to compute the derivative of V(t) along the trajectories of the 
prediction error dynamics (10) and the adaptive laws (9). Using the properties of the projection operator 
(see for example Ref. 8 ) we obtain the inequality 

V(t) < -2A y 2 (t) + 2 y(t)5 2 {t) - 6 T (t)QS(t) . (12) 

Completing the squares we arive at 

V(t) < — (2A - 1 )y 2 {t) - (A min (Q) - l)||<5(t)|| 2 . (13) 

Choosing A m i n (Q) > 1 and 2A > 1, we conclude from the LaSalle- Yoshizawa theorem (, 6 p.24) that y(t), 
d m (t) and 1 9(t) are globally uniformly bounded, and y(t) — >■ 0 as t — > 00. In particular, there exists /?i > 0 
such that \dm(t)\ + II^WII 2 < /3 2 - 

The next lemma gives the bound on the output prediction error. 
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Lemma III. 2 y(t) satisfies the bound 


\m\<p2e- Vot 


fa_ 


where fa = yJ\2V{Q) — ^-| and vo = min(A m i n (Q) — 1, 2A — 1). 

Proof. It is easy to see that the inequality (13) can be written as 

V(t) < -uoV(t) + ^/3 1 , 

7 


which implies that 


V(t)< 


V(0) 



ft 

2 7 ' 


Recalling that ||y(t)|| 2 < 2 V(t), we obtain 


lls/(*)ll < 





(14) 


(15) 


(16) 


(17) 


The bound (14) follows from the fact that \J a + b < faa + y/b for any a > 0, b > 0. 

When y{t) and u(t) are bounded (which will be provided by the control design), a tighter bound on y(t) 
can be derived through the bound on the adaptive error signal f]{t), which is given by the following lemma. 

Lemma III. 3 Let the estimates y(t ) and i9(f) be generated by the identification model (8) and (9). In 
addition, let y(t) and u(t ) be bounded. Then fj(t) and y(t) satisfy the following bounds 

\v(t)\ < fae~ vt + ^ 

\y(t)\ < fae~ vt + ^ , 

where the constants fii > 0, i = 3, . . . , 6 and v i > 0 are defined in the proof. 

Proof. Following Ref., 10 it is straightforward to show that fj(t ) satisfies the dynamics 

fjfit) + A fji(t) + 'yp{t)fn(t ) = -7p(f)y(f) + cr a (t ) + f a (t) - 7 p{t)S 2 (t) , (20) 

where p(t) = wf, 2 {t) + h T (t)h(t), r a {t) = & m (t)u> mi2 (i) + i) T (t)h(t). Since y(t) and u{t) are bounded, and 
f(y), $(j/) are smooth functions, s(t), fat) and fat) are bounded. These imply that w m fat), h(t) and y(t) 
are bounded. Also, it follows from (6) that y{t) is bounded, which implies that Wmfat), fat ) are bounded 
as well. Therefore, there exist positive constants c\, C2, C3 such that ||p(t)||,c 00 < Ci, || / o(t)||,c oo < C2 and 
ll r a(^)IUco < C3. On the other hand S 2 (t) satisfies the bound |<52 (^) | < C 4 e ~ Ult , where C4 > 0 depends on the 
initial conditions and Z7 is the decay rate of the Hurwitz matrix A 0 . Then, it follows from the results of 
Ref. 10 that choosing A > 2facfiy damps the oscillations in 77(f) and guarantees the bound 

\m\<foe- vt +c 2 \y{f)\ + ^\r a (t)\, (21) 

where v = min(zT), v-l, ^2), ^2 is proportional to fay , and the positive constants fa and C4 are independent 
of 7 (see details in 10 ). Substituting (14) we arrive to (18) with fa = c 2 fa + C3C5. 

A tighter bound on y(t) is obtained from (10) by direct integration. 

\y(t)\ < yh i e ~ Vt - e_A ‘] + ffaTfa l e ~^ - e_A ‘] + -Q= [! - < fae~ ut + , ( 22 ) 

since > v by definition and A > v for fast adaptation. 


(18) 

(19) 
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IV. Control Design 


In this section, we first design a conventional backstepping controller for the known system (ideal con- 
troller), which is used only for analysis purposes, then design the command filtered certainty equivalence 
controller, which is actually being implemented. 


A. Known System 

We follow the standard backstepping procedure from Ref. 6 and design the ideal controller for the system 
(7) ignoring the term S 2 {t). The reason is that the actual control signal is defined using the identification 


model, which does not evolve that term. Let the error variables be defined as 

e\{t) = y\t)-y r {t) (23) 

= lWi * = 2, •••,?, (24) 

where the stabilizing functions o®(t), * = 1, . . . , q — 1 have the form 

ot°{t) = a^[-ce?(i) - fi(y°) - s%(t) - h 0T + x r ,2 (*)] (25) 

a °(t) = ~d m e\{t) - c 2 e° 2 (t) + k 2 W° m l (t) + d?(f) (26) 


“i (t) = -e°_i(i) - Cie°(t) + kiW° ml (t) + a^t), i = 3,...,q-l. 


with Ci > 0, i = 1, . . . ,q— 1 being design parameters. Here and in the following derivations, the superscript 
0 indicates that the variables correspond to the ideal case. We define the ideal control signal as 


u°(t) = -e9_!(t) - c q e°(t) - w^ q+1 (t) + k q w° m ^{t) + d°_ x (t) , 
It is easy to see that the error system has the form 

e°(i) = A id e°(t ) , 


(27) 

(28) 


where 


Aid = 


-Cl d m 0 0 

dm ^2 1 0 

0 -1 -c 3 1 


0 

0 

0 


.0 0 0 0 ... —c q 

Therefore it is exponentially stable. That is, y°(t) = y r (t ) + e 5 6 (t) exponentially converges to y r (t). In 
addition, from r{t) £ Coo it follows that x r (t) £ Coo-, which along with e°(t) £ Coo can be used to recursively 
show that a 6 (f), d 6 (f), * = 1, . . . , q — 1, w m ,i, i = 1, . . . , q + 1 are bounded. Then, u°(t) £ Coo is bounded, 
and the boundedness of and hence all closed-loop signals follow. 


B. Command Filtering Certainty Equivalent Control 

Here, we use command filtering approach from Ref. 2 and design a certainty equivalence control for the 


identification model augmented with the last q — 1 equations in (7). 

The uncompensated errors are introduced as 

ei {t) = y{t) - y r (t) (29) 

ej(f) = w mii (t) - 0 -j-i,i(t), i = 2,...,q, (30) 

where <7j_i.i(i) is the first state of the second order LTI system, which is used to filter the stabilizing 
functions. Here, these functions are defined as 

ai (t) = y^— [— cei (f ) - /i (y) - s 2 (t) - h T (t)tf + x r ,2 (i)] (31) 

dm (t) 

o 2 (t) = -d m (t)ei(t) — c 2 e 2 [t) + k 2 w m ,i(t) +uja lt2 (t) (32) 

Oi{t) = — ei_i(t) - Cidit) + hwmAt) + uj(Ti—i 2 (t) , i = 3,...,q-l, 
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where cq-i^t) is the second state of the filter’s dynamics, which is given by the equations 

= cja i}2 (t) (33) 

= -2Cwc - w[c T itl (t) - Oti(t)} 
i = 1, . . . , n — 1 , 

with the initial conditions CTi.i(O) = aj(0) and 2 ( 0 ) = 0. 

The compensated errors are introduced as e c (t) = e(t) — £(t), where £(f) is dynamically defined as 

ii{t) = dij(t) — a\(t) + (34) 

£i(t) = &i, i(t) - oti(t) + £i+i{t), i = 2,...,n-l 

with £;(0) = 0 for * = 1, . . . , n — 1, and £ n (t) = 0. These modifications lead to the definition of the actual 
control signal as 

u(t) = -e\_ i(t) - c q e q (t) - w m , 9+ i (t) + k q w m>1 (t) + u(T q _ ia (t) . (35) 

The resulting compensated error dynamics are given by the equation 

e c [t) = A{t)e c {t) + A c q y(t) , (36) 

where A(t) is obtained from A ^ by replacing d m with its estimate d m (t). Obviously, the error system 
(36) is asymptotically stable, since y(t) — ► 0 as t — > 00 according to Section III. On the other hand, the 
uncompensated error e(t) has dynamics 

e(t) = A(t)e(t ) — a(t) + &(t) + \c q y(t) . (37) 


where 



dm (t)ai (t) 



a(t) = 

a q -i(t) 

, a(t) = 

Vq-1,1 (t) 


0 


0 


The performance of the designed controller is given by the following theorem, the proof of which can be 
carried out following the steps from Ref. 11 

Theorem IV. 1 Let the controller be defined according to identification model (8), adaptive law (9), and 
command filtered backstepping scheme given by (29), (31), (33), (34) and (35). Then the input and output 
tracking errors satisfy the following upper bounds 


\u(t)-u 0 (t)\<fi 7 e-^ + ^+O{e) (38) 

I y(t) - Vr(t) I < /3 9 e- ! ' 3t + + O (e) , (39) 


where (3 7 , fis, fig, fiio are positive constants independent of 7 , vg is the decay rate of the reference model, 
e = 1 /u) (the proper choice of ( and ui is discussed in Ref. 2 ), and the notation O (e) is adopted from Reffi 
(p. 383). 


V. Aerospace Applications 


A. Control of Wing Rock 

We consider a wing rock motion of a slender delta wing, which is given by the equation 


4>{t) = /(</>, 4>) + b a S a (t ) . 


where (j>(t) is the roll angle, 5 a (t) is the aileron deflection and f(<f>,(j>) has the form (see Ref. 9 for details) 


f{fi,<fi) = Oxfi + + 9^\(t) 2 \(j) + 04\fi\(j) + = 6 T ip((j),p) . 


(40) 
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Modeling ailerons as a first order dynamics 

rS a (t) = u(t)-d a (t), (41) 

where r is the ailerons time constant, and introducing a new variable x 3 (t) = b a S a (t ), we represent the 
system in the state space form as 


±i(t)=x 2 {t) (42) 

±2 (t) = x 3 {t) + 0 T ip(x i, x 2 ) 
x 3 (t) = — kx 3 (t ) + bu(t ) , 

where x\{t) = <j>(t) denotes the roll angle, x 2 (t) = p[t) denotes the roll rate, k = \ is assumed to be known, 
b = ^ is assumed to be unknown, and u(t) is the control signal to be designed assuming that Xi(t) and 
x 2 (t) are available for feedback. 


Systems response to ideal command filtered backstepping control 



time 


Figure 1. Step response to the ideal command filtered output feedback control. 


The system consists of two blocks. First contains only the first equation, which is free of uncertainties, 
the second consists of the last two equations, which contain unknown parameters. Since x 3 is not accessible, 
we transform the second block into the parametric output feedback form using a linear transformation 

z 2 (t) = x 2 (t) 

z 3 (t ) = kx 2 (t) + x 3 (t) . (43) 

The transformed system has a form 

zi(t) = z 2 {t) (44) 

z 2 {t) = z 3 (t) - kz 2 (t ) + 9 T ip(x i,x 2 ) 
z 3 (t) = bu(t ) + kd T cp(xi,x 2 ) , 

where for uniformity we also introduce Zi(t) = X\ (t). Using the filtered transformation (46) we write the 
system as 

zi(t) = z 2 (t) (45) 

z 2 (t) = b/j, 2 (t) + s 2 (t) - kz 2 (t) + i9 T h(f) 
fi 2 (t) = u(t) - k 02 fxx{t ) , 
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Systems response to cerainty equivalent control 




Figure 2. Step response of the certainty equivalent control. 


Aileron deflection 



Figure 3. Actual aileron deflection and certainty equivalent adaptive control signal vs ideal ones for the step 
command. 


where h(t) = [S( 2 )(i) 4- p> T (x i,x 2 ) ' 02 (^)] T i and s(t), 3(i), fi(t) are generated according to equations 


s(t) = A 0 s(t ) + 

1 

o 

i 

x 2 (t), E(t)A 0 E(t) + 

ip T {x i,X 2 ) 

7 1 / \ 

, ii(t) = A 0 fi(t) + 

0 


&0,2 


kip (X 1 ,X2) 


1 


u(t), (46) 
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where 


An = 


— ko,i 1 
^- 0,2 0 

The identification model is introduced for the ^-dynamics, where the uncertainties are lumped in. 


Z2 (t) = b(t)y, 2 (t) + s 2 (t) - kz 2 {t) + 1? ( t)h(t ) + A [z 2 {t) - z 2 (t)] 
where the estimates b(t) and $ (t) are generated according to adaptive laws (9). 

Systems response to ideal command filtered backstepping control 



(47) 



Figure 4. Sinusoidal response of the ideal command filtered output feedback control. 


In the first step of the control design we introduce the error ei (t) = y(t) — y r (t) = Z\{t) — x r ,i{t) and the 
stabilizing function a\(t) = —Ci{t)e\(t) + x r ^ 2 , which results in 


ei(t) = -ciei (t) + e 2 (t) + z 2 (t ) , 


(48) 


where e 2 (t) = z 2 (t) — a\(t). Since differentiation of a\(t) involves only linear terms, we do not use command 
filtering in this step. This results into the equation 


e 2 (t) = b(t)n 2 {t) + s 2 (t) - kz 2 (t ) + ( t)h(t ) - x ry3 - + cie 2 (i) + (A + Ci)z 2 (t) . 


(49) 


The stabilizing function a 2 (t) is modified to include a term — e 3 (t) to counteract the effect of the error term 
e 2 (t) in (48), and has the form 

a 2 (t) = -J— (cf - l)ei(t) - (ci + c 2 )e 2 (t) - s 2 (t ) + kz 2 (t) - ^ {t)h(t) + x Vy3 
b (t ) L 

The last error signal is defined as e 3 (t) = y, 2 (t) — o\ (t), where a i (f) is generated by the system 

&l (t) = wcr 2 (t) (50) 

<j 2 (t) = -2(uj(j 2 (t) - w[ai(t) - a 2 {t)] . 

This reduces the error dynamics (49) into 

e 2 (t) = -ei (t) - c 2 e 2 (t) + b(t)e 3 (t) + b(t ) [a ± (t) - a 2 {t)] + (A + ci)z 2 (t) . (51) 


9 of 14 


American Institute of Aeronautics and Astronautics 


Systems response to cerainty equivalent control 




Figure 5. Sinusoidal response to the certainty equivalent control. 


Aileron deflection 



Aileron deflection command 



Figure 6. Actual aileron deflection and certainty equivalent adaptive control signal vs ideal ones for the 
sinusoidal command. 


The control signal is defined as 

u(t) = -b(t)e 2 (t) - c 3 e 3 (t) + k 02 /j.i(t) - ujcr 2 (t ), 
which reduces the dynamics of e 3 (t) into 

e 3 {t) = — S(t)ea(t) - c 3 e 3 (f) (52) 

We test the algorithm in simulations for the wing rock model from Ref. 9 with 6\ = —0.0186, d 2 = 0.0152, 
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0 3 = —0.0625, 04 = 0.0095, d 5 = 0.0215, b a = 1 and r = 1/15. We choose k 0 = [3 2] T for the filtered 
input and output transformation. For the reference model we set k = [4.86 7.56 4.38] T and k r = 4.86]. The 
command filtering is implemented with u> = 500 and £ = 0.8. We run simulations from the initial conditions 
</(0) = 6 deg, p( 0) = Meg /sec, 5 n (0) = 0 deg. The adaptation rate is set to 7 = 25000 with A = 2^/7, and 
the control gains are set to Ci = C2 = C3 = 2. The lower bound for b(t) is set to 0.01. 

Figure 1 displays the tracking performance of the ideal command filtered output feedback backstepping 
control for the bank angle step command of magnitude 15 deg, assuming that the parameters are known to 
the controller. The performance of the certainty equivalence control with respect to the ideal response is 
presented in Figure 2, and the corresponding aileron deflections and commands to the aileron actuators are 
presented in Figure 3. It can be seen that the presented output feedback M-MRAC achieves good tracking 
performance for both output and control signals. 

Next set of figures (Figures 4,5 and 6) represent the systems performance for the sinusoidal external 
command of magnitude 15 deg and frequency of 1 rad/ sec. Figure 4 displays the tracking performance of the 
ideal controller with respect to the reference model, and Figure 5 displays the performance of the certainty 
equivalent control with respect to ideal one. Again, it can be observed that the M-MRAC exhibits adequate 
behavior as in output tracking, as well as in the input tracking which is apparent from Figure 6. 


B. Control of Short Period Dynamics 


As an example we consider the longitudinal control of a damaged transport class model (TCM) as defined in 
Ref. 3 The short period dynamics of the GTM can be represented as (see for example Ref. 12 for definitions) 

a(t) = ^ a(t ) + ^1 + ^ q(t ) (53) 

q(t) = M a a(t) + M q q(t) + M s J e {t) , 

where a is the angle of attack, q is the pitch rate and S e is the elevator deflection, the simplified dynamics 
of which are modeled as a first order system 

5 e (t) = k e [u(t) - 8 e (t)} , (54) 


where u(t) is the control input to the elevator actuator, k e is the actuator’s gain. We assume that the 
controlled output is the pitch rate, which is the only measurement available for feedback. 

Deriving the transfer function from u to q we obtain 


G U q[s) — 


1 - fcey 


S 3 + (ke - Mq - S 2 + 


Mq^-M c 


(l + 1 / ) — (M q + 


Mq^-M c 


{}+¥) 


therefore the system can be easily transformed into the observer form 


±i(t) = x 2 (t) + diXxit) (55) 

±2 (t) = x 3 (t) + d 2 Xi (t) + biu(t) 

X 3 (f) = 0 3 xi{t) + b 0 u(t) 

y[t) = x 1 {t), (56) 


where the form of the unknown parameters 9\ = k e — M q — ^ , 9 2 = M q ^y- — M a ^1 + — k e (M q + Z y ) , 

03 = k e Mq^- — M a ^1 + , b\ = k e and b 0 = -k^-. Following Ref. 6 (p. 421) we use the input and 

output filters for the special case of linear systems 


S u (t) = A 0 s u (t) + 

l 

O O 

1 

u(t), Sy(t) = A 0 s u (t) + 

1 

0 0 
1 

II 

0 

'IS 

— fco.i 1 0 

— fco,2 0 1 


1 


1 


1 

O 

O 

CO 

0 

1 

1 


(57) 


to generate vectors Wj(t) = A 3 Q s u (t), j = 0,1, s(t) = AQS v (t) and £-(t) = A J 0 s v (t), j = 0,1,2 and the 
matrix S(f) = [£ 2 (£) £1 (i) £ 0 (*)]- Then the system (55) can be written as 

y{t) = 6iwi, 2 (0 + s 2 (t) + h T (t)'d (58) 

wi ,2(t) = u(t) + wi i3 (t) - k 0 , 2 Wl,l(t) , 
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Systems response to ideal command filtered backstepping control 



Ideal command filtered backstepping control 



Figure 7. Ideal pitch rate response and ideal command filtered output feedback control. 


Systems response to cerainty equivalent control 




Figure 8. Pitch rate response and the corresponding angle of attack. 


where h T (t) = [wo^(t) S( 2 ) + y{t)[ 1 0 0]], and = [&o #i 62 ^] T - The system is now in the from (7), 
therefore the presented output feedback M-MRAC can be applied. 

For the simulation the numeric values of matrices A and B are computed using the vortex-lattice method 
(see Ref. 7 for details), for the straight level flight condition at an altitude of 30000/t and a cruise speed of 
0.8 M 

The damage scenario corresponds to 60% loss of the left horizontal tail, which roughly corresponds to 59% 
loss of pitch control effectiveness, 67% loss of pitch damping and 59% loss of pitch stiffness. It is assumed 
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Elevator deflection 



Figure 9. 


Actual elevator deflection and certainty equivalent adaptive control signal vs ideal ones. 


that the aircraft was in steady level flight condition at an altitude of 30000/f and a cruise speed of 0.8 M 
before the failure. 

We choose fco = [6 11 6] T to generate AO, a unity gain second order reference model with frequency of 
2.5 and damping ratio of 0.8 and implement the command filtering with a filter, which has ui = 300 and 
£ = 0.8. We run simulations from the initial conditions q( 0) = Odeg/sec, a(0) = 3 deg, <5 e (0) = Odeg. The 
adaptation rate is set to 7 = 15000 with A = 2^7, and the control gains are set to c\ = C2 = 0.5. The 
upper bound for 61 (t) is set to —0.01 (b\ is known to be negative). The external command is a square wave 
of magnitude 5 deg /sec and of frequency 7r/8 rad/ sec. In order to make the command differentiable, it is 
filtered trough a first order filter 

Figure 7 represents the tracking performance and the ideal command filtered output feedback backstep- 
ping control signal. The first plot of Figure 8 displays the output tracking performance of the certainty 
equivalence M-MRAC with respect to the ideal response. It can be observed that M-MRAC provides a good 
output tracking. The second plot in the same figure displays the time history of the angle of attack. It 
can be seen that it closely follows the one generated trough the ideal control. The corresponding elevator 
deflection and the control input to the actuator are presented in Figure 9. It is seen that the output feedback 
M-MRAC achieves good tracking for the control signal as well. 

VI. Concluding remarks 

We have presented an adaptive output feedback backstepping control method for nonlinear systems with 
unmatched uncertainties that follows the certainty equivalence principle without over-parametrization. The 
approach uses a fast identification model, which is independent of the control design. This separation of 
the parameter estimation and control design enables the designer to achieves desired transient and steady 
state properties. In the control design process we utilize the command filtered backstepping approach. The 
presented method was applied to two aerospace control problems with good results. 
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